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A  CODE  FOR  THE  SECONDARY  ELECTRON  ENERGY  DISTRIBUTION 
IN  AIR  AND  SOME  APPLICATIONS 

1.  INTRODUCTION 

The  interaction  of  high  energy  electrons  (>1  KeV)  with  atoms  and 
molecules,  the  generation  of  the  secondary  electrons  and  their  subsequent 
energy  degradation  is  a  subject  of  great  interest  in  many  areas  of  applied 
physics.  Among  these  areas  are:  precipitation  of  energetic  electrons  in  the 
upper  atmosphere  resulting  in  auroral  phenomena,  production  of  laser  radiation 
by  high  energy-high  current  electron  beams,  and  electron  beam  generated 
discharges. 

Various  methods  have  been  utilized  to  address  the  electron  energy 
deposition  problem.  These  may  be  classified  in  different  ways  depending  on 
one's  emphasis.  Major  classifications  are  as  follows:  1)  Transport  versus 
local  descriptions,  2)  Continuous  slowing  down  approximation  (CSDA)  versus 
discrete  energy  loss  description  (DEL) ,  and  3)  Transport  or  local  energy  loss 
equation  versus  Monte  Carlo  approach,  there  is  an  extensive  literature  on 
this  subject  as  it  relates  to  applications  cited  above.  Here  we  will  just 
note  some  of  the  works  representative  of  past  efforts.  Green  and  his 
colleagues  '  have  approached  the  problem  through  the  continuous  slowing 
down  approximation  and  have  applied  their  calculations  to  auroral  phenomena  as 
well  as  deposition  in  other  gaseous  media.  The  continuous  slowing  down 
approximation  has  been  improved  to  include  discrete  energy  loss  for  electrons 

4  5 

with  energies  below  SOOeV  and  has  been  applied  to  auroral  phenomena  '  as 

6*8  9  10 

well  as  other  gaseous  elements  .  Monte  Carlo  and  Fokker-Planck  methods 
have  also  been  utilized  to  calculate  the  energy  deposition  as  well  as  the 
spatial  distribution  of  primary  elctrons.  The  steady  state  forms  of  the 
Boltzmann  equation  provide  another  method1  1 2 , 1 3  to  calculate  the  equilibrium 
distribution  of  the  secondary  electrons. 

Manuacript  approved  September  30, 1982. 
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In  this  work,  we  are  specifically  interested  in  detailed  energy 
partitioning  in  air  during  and  following  the  passage  of  a  high  energy  electron 
beam  through  the  air.  For  this  purpose  we  have  developed  Code  SED  which 
provides  a  discrete  energy  loss  description  for  specification  of  the  secondary 
electron  distribution  function.  The  approach  is  basically  that  given  in 
reference  11  but  limited  to  local  deposition.  This  code  has  been  coupled  to  a 
time  dependent  air  chemistry  code14  (CHMAIR-II)  which  describes  the  many 
effects  on  neutral  and  ionized  species  during  and  following  electron  beam 
deposition.  In  this  report,  however,  we  limit  our  discussions  to  the 
secondary  electron  distribution  obtained  using  code  SED. 


2.  APPROACH 

We  wish  to  specify  the  quantity  ne(E,t)  which  is  the  density  of  secondary 
electrons  per  unit  volume  per  unit  energy  (cm~3eV“1).  Given  ne,  we  may  obtain 
various  production  rates  of  interest  given  by  P: 

/•  -3  .1 

P(t)  »  n/  ne(E,t)  v(E)  a  (E)  dE  cm  sec  (1) 

where  n  is  the  density  of  the  target  species  and  o  is  the  electron  impact 
cross  section.  Most  of  the  contributions  to  P  comes  from  below  100  eV  and  for 
this  reason  we  restrict  ourselves  to  low  energies  for  which  transport  effects 


are  unimportant 


The  density  nfi  is  specified  using 


OM  i 

-  S(B,t)  -  nlal  <E)v(E)ne(E,t) 

+n  (t)  r|—  [L  (E)v(E)  n  (E,t)l 
p  3E  1  p  e  1 


4k 


(E',E)v(E')n  (S'* ,t)dE‘* 
e 


(2) 


where  S  is  the  volume  production  rate  in  e/cm3-s-eV,  n^  is  the  1th  species 
causing  energy  loss,  the  o's  are  energy  loss  cross  sections,  is  the  plasma 
density,  and  Lp  is  the  loss  function  for  loss  to  the  plasma  electrons.  The 
species  treated  in  this  equation  are  Nj  and  02*  The  electrons  lose  but  do  not 
gain  energy  in  the  description,  thus  allowing  us  to  obtain  ne  one  energy  at  a 
time  starting  at  the  high  energy  end.  The  energy  dependence  of  nQ  within  the 
integral  is  specified  in  a  manner  similar  to  that  used  by  Strickland,  et  al 1 1 . 


3.  SPECIFICATION  OF  S(E,t) 

The  secondary  electron  volume  production  rate  S(E,t)  is  given  by 

S(E,t)  -  W  Snjt<Je(Eb,E)  (3> 

l 

-2  - 1 

where  is  the  beam  current  density  in  e  •  cm  •  s  and  o^(E^,E)  is  the 
differential  ionization  cross  section  for  species  l  .  We  use  the  expression 
by  Porter,  et  al3  for  this  cross  section  which  is  similar  to  the  Moeller15 
formula  for  large  momentum  transfer  values.  Figure  1  shows  the  differential 
ionization  cross  section  of  Nj  from  1  to  lOOOeV  for  Etj  *  1  MeV.  The  Moeller15 
cross  section  is  included  to  conpare  with  the  Porter  expression  at  large  E_ 
values  where  both  are  in  close  agreement. 


ELECTRON  IMPACT  CROSS  SECTIONS 


In  the  discrete  energy  loss  description,  cross  sections  are  required  for 
all  important  inelastic  scatterings.  This  amounts  to  approximately  ten  cross 
sections  per  species,  most  of  which  apply  to  the  process  of  excitation. 
Ionization  accounts  for  the  remaining  cross  sections.  For 
excitation,  c^(E",E)  may  be  expressed  as 

<J£k(E',E>  -  «[E'  -  (E+W^)]  (4) 

where  the  delta  function  restricts  energy  loss  to  the  excitation  energy 
W  £k  .  For  ionization,  we  use  the  form 

0B.  (E',E)  -  a.A  E')  p(E',E)  (5) 

where  E")  is  the  total  ionization  cross  section  in  cm  for  ion  state  k 

and  p  is  a  function  with  normalization 

-Wk)/2 

E“*  ,E )  dE  -  1  •  (6) 

o 

The  integration  is  over  the  secondary  electron  energy  defined  to  be  the  lesser 
of  the  energies  of  the  two  electrons  emerging  from  an  ionization  event.  The 
form  of  p  is  chosen  to  be  that  used  by  Strickland,  et  al11. 

The  individual  cross  sections  will  not  be  presented  here.  They  can  be 
found,  however,  in  a  recent  conqpilation  by  Ali16.  The  total  cross  section  for 
each  species  may  be  seen,  however,  in  Figure  2.  The  N2  feature  peaking  at  2.5 
eV  results  from  vibrational  excitation  of  the  ground  state.  This  single 
feature  is  actually  composed16  of  eight  cross  sections  for  excitation  through 
vibrational  level  v"  ■  8.  The  minimum  in  the  N,  cross  section  near  6  eV  will 


be  seen  to  have  a  noticeable  effect  on  the  secondary  electron  distribution 
near  this  energy.  The  distribution  takes  on  large  values  in  this  region  since 
electrons  find  themselves  trapped  in  velocity  space  for  a  time  due  to  the 
small  inelastic  collision  frequency. 


5.  ENERGY  LOSS  TO  THE  PLASMA  ELECTRONS 

Equation  (2)  includes  an  expression  for  energy  loss  to  the  plasma 
electrons.  This  term  becomes  important  for  low  energy  secondaries  especially 
when  fewer  electronic  state  excitation  processes  occur  with  decreasing 
energy.  For  L  ,  we  are  using  the  following  expressions  from  Schunk  and  Hays1^ 


L  (E) 
P 


1.3  x  IQ"1 3 
E 


8.2  x  109  E1.5 

n  z 
P 


E  <  20  eV  (7) 


1.3  x  10~13 
E 


In 


5.4  x  1010 


n* 


E  ,  E  >  20  eV 


P 


(8) 


The  plasma  electron  density  np  appears  explicitly  in  the  plasma  loss  term  in 
equation  (2)  and  in  Lp  above.  For  our  time  dependent  results,  np  is 
calculated  self-consistently  along  with  n@.  For  our  steady  state  results  to 
be  presented  below,  np  is  simply  an  input  parameter  to  the  calculations. 


6.  MATRIX  APPROXIMATION  TO  THE  SECONDARY  ELECTRON  EQUATION 
Equation  (2)  can  be  written  as 


(9) 


where  And  l ^  are  the  volume  production  and  loss  rates  for  species  i.  In 
the  case  of  secondary  electrons,  the  index  i  refers  to  energy  E.  To  arrive  at 
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Equation  (9),  approximations  oust  be  made  to  the  integral  and  plasma  loss 
terms  in  Equation  (2) .  This  is  done  by  first  specifying  an  energy  grid  over 
the  range  of  interest.  The  density  n@  within  the  integral  is  then  given  a 
quadratic  dependence  over  a  given  energy  zone  following  the  method  of 
Strickland,  et  al11.  This  leads  to  replacement  of  the  integral  by  a  sum  of 
terms  of  the  form 


i 

E 

Z  j-1 


R.  .  n  . 
ij  ej 


For  the  plasma  term,  one  sided,  finite  differencing  is  used  over  the  interval 
(E^,  E^^)  where  Ei_1  is  greater  than  Ei#  Equation  (2)  may  then  be  expressed 
in  the  form 


3n 


ei 


at 


l  ii 


+  n 


-Ei 

AE 


i  \ 

-  v.  I  n  . 

i  L)  61 


(10) 


i-1  J l 


L  ,  .v. 


:  _  P, i-i  i-i 

+  Z  n,  I  ,  R.  .  n  ,  +  n  -  n  .  . 

K  l  3-1  ij  ej  p  AE^  e,i-1 


+  s. 


The  matrix  element  R^  is  obtained  from  combinations  of  terms  whose  integral 


parts  have  the  form 


E 


j-1 

J 

j 


J*  alk  (E",E)  v(E")  f(E")  dE" 


where  f(E  ")  is  either  In^E  ",  In  E  ",  or  1. 

7.  ENERGY  CONSERVATION 

Two  tests  have  been  applied  to  our  solution  of  Equation  (10)  in  steady 

state.  The  first  determines  how  well  n  (t)  conserves  energy  and  in  the 

© 


second,  we  calculate  the  energy  expended  to  generate  an  ion  pair,  A  W. 


Th 


accepted  value  of  A  W.  for  air  Is  ~  35  eV/ion  pair,  'Hie  conservation  test 

*P 

determines  how  well  the  following  equality  holds: 


E  W 
f,k 


Ik 


ik 


%/' 


L  (E) 
P 


v  n 


(E) 


S(E)  E  dE 


(11) 


The  conservation  principle  states  that  the  rate  of  energy  transferred  to  the 

ion  and  excited  states  and  to  the  plasma  equals  the  rate  of  energy  going  into 

the  secondary  electron  spectrum.  TSie  production  terms  P  refer  only  to 

production  by  secondary  electrons.  We  find  that  the  solution  n  (E)  leads  to 

© 

less  than  10%  error  in  satisfying  Eq.  (11). 

For  AW.  ,  we  are  obtaining  a  value  of  between  34  and  36  which  agrees 

ip 

well  with  the  accepted  value.  To  calculate  AW.  ,  we  first  determine  how 

xp 

much  energy  is  deposited  in  a  unit  volume  per  unit  time  by  passage  of  the 
primary  electrons  through  it.  The  Porter3  cross  section  is  used  to  give  the 
total  primary  ionization  rate  which  in  turn  gives  the  energy  going  into 
potential  energy  of  the  ion  states.  This  same  cross  section  gives  S,  which 
when  multiplied  by  E  and  then  integrated  gives  the  rest  of  the  beam  energy 
loss  by  ionization.  There  is  also  energy  loss  to  excited  states.  Based  on 
cross  sections  in  the  keV  regime,  the  number  of  excitation  events  will  be 
similar  to  the  number  of  ionization  events  caused  by  direct  beam  electron 
impact.  This  provides  us  with  the  loss  due  to  excitation  once  given  the 
average  excitation  energy  threshold. 

The  other  quantity  needed  to  specify  A  W^  is  the  total  number  of  ion 
pairs  produced  by  the  beam  energy  loss.  The  integral  of  S(E)  gives  the  number 
produced  by  direct  beam  ionization.  Equation  (2)  gives  the  number  of  ion 
pairs  by  all  generations  of  secondary  electrons  for  each  ionization  process. 
Summing  over  processes  and  adding  the  result  to  the  direct  beam  contribution 
gives  the  desired  result. 
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8.  RESULTS  AND  DISCUSSION 

The  problem  being  addressed  is  the  incremental  energy  loss  of  an  electron 
beam  on  its  passage  through  some  thickness  6s  .  We  wish  to  know  how  the 
deposited  energy  is  distributed  among  various  ion  and  excited  states.  The 
production  rates  for  these  states  may  be  divided  into  two  parts  -  one  produced 
by  direct  impact  of  the  beam  electrons  and  the  other  by  inpact  of  secondary 
electrons.  This  requires  the  inpact  cross  sections  both  at  the  primary  energy 
and  over  che  energy  range  important  to  secondary  electron  inpact  processes. 

Calculations  have  been  performed  for  energy  deposition  in  air  with 
primary  electron  fluxes  having  energies  of  1,  10,  100,  and  1000  keV.  The 
easiest  part  of  these  calculations  is  the  specification  of  P^,  the  production 
rates  by  primary  or  beam  electron  inpact.  Pfa  is  given  by 


The  production  rate  produced  by  secondary  electron  inpact  Pg  was  given 
previously  by  Equation  (1).  The  needed  secondary  electron  flux  vn  has  been 
obtained  from  Equation  (2)  in  its  steady  state  form  with  the  source  spectrum 
given  by  Equation  (3).  The  time  independent  results  are  satisfactory  here 
since  our  aim  is  to  show  how  energy  becomes  partitioned  among  the  various 
states.  For  chemistry  modeling  purposes,  however,  explicit  time  dependence  is 
necessary. 

Before  presenting  results  from  solutions  of  Equation  (2) ,  we  wish  to  note 
how  the  energy  is  distributed  in  the  secondary  electron  source  spectrum.  We 
are  thus  discussing  the  behavior  of  the  applied  differential  ionization  cross 
sections  for  N2  and  02«  The  purpose  for  presenting  this  information  is  to 
note  how  much  source  energy  is  not  taken  into  account  when  Equation  (2)  is 
solved  over  less  than  the  fully  allowed  secondary  electron  energy  range.  This 
is  pertinent  to  our  calculations  when  the  primary  energy  is  above  a  few 


10 


hundred  keV  since  we  do  not  account  for  secondaries  with  energies  greater  than 

100  keV.  Figure  3  shows  the  fraction  of  the  initially  generated  secondary 

electron  energy  f_(E_)  at  energies  less  than  E_  for  primary  energies  of  1,  10, 

100  and  1000  keV.  In  each  case,  the  cutoff  in  f_  comes  where  E_  is  half  of 

s  s 

the  primary  energy.  This  defines  the  upper  limit  of  Es.  We  observe  that  the 
high  energy  tail  of  the  secondary  electron  source  spectrum  does  contain 
measurable  energy  in  spite  of  the  fact  that  the  electrons  available  are 
insignificant  conqpared  to  the  total  number  within  the  spectrum.  Figure  4 
recasts  the  information  from  Figure  3  in  terms  of  the  average  energy  <E  >  of 
the  initially  produced  secondary  electrons.  At  e  given  energy  E_,  <E„>  is  the 

9  9 

average  from  O  eV  to  E  .  The  rise  in  <E  >  with  increasing  energy  is  somewhat 

3  9 

steeper  than  would  be  obtained  using  the  Moeller  expression  for 

-2 

c(E,  ,E  )  which  varies  as  E  .  This  is  illustrated  in  Figure  4  by  means  of 
D  S  S 

the  dashed  curve  which  shows  how  <Eg>  varies  for  E^j  ■  1  MeV  assuming  the 
-2 

Eg  dependence.  We  now  proceed  to  discuss  our  results  but  shall  refer  back 
to  Figures  3  and  4  for  those  calculations  which  have  not  taken  into  account 
all  available  secondary  electron  energy. 

Examples  of  the  calculated  flux  vne(E)  are  shown  in  Figure  5  for  E^  *  1 
Mev.  For  these  results,  the  maximum  energy  treated  was  20  keV.  Referring 
back  to  Figure  3,  we  are  thus  accounting  for  60%  of  the  available  secondary 
electron  energy  based  on  the  use  of  the  Porter  et  al.  ionization  cross 
section.  The  calculated  flux  has  been  divided  by  the  primary  flux  thus  giving 
the  units  of  e/eV-Sp.  The  spectra  shown  are  for  pure  Nj  and  pure  O2,  in  each 
case  at  standard  density.  The  spectrum  for  air  is  similar  to  that  of  Nj.  The 
point  of  calculating  the  exhibited  spectra  was  to  determine  the  sensitivity  of 
the  secondary  spectrum  to  cooposition  changes.  This  is  of  interest  to  us  in 
our  chemistry  modeling  since  we  do  observe  compositional  changes  in  the  major 


Fig.  3  —  Fraction  of  secondary  electron  source  energy,  fg  between  O  and  E,  for  primary 
energies  as  shown,  f,  is  based  on  the  ionization  cross  section  of  Porter  et  al3 . 


Fig.  4  —  Average  initial  secondary  electron  energy  <ES)  obtained  using 
<E,)  -  <E,>max  f,(E,)  where  is  the  average  energy  over  the  full 

range  of  secondary  electron  energies.  The  solid  curves  were  obtained 
using  the  ionization  cross  section  of  Porter  et  al3  while  the  dashed 
curve  was  obtained  using  the  Moeller  cross  section. 


E,  (eV) 

Fig.  5  —  Secondary  electron  fluxes  vn^E)  for  pure  N2  and  pure  02  compositions 
at  standard  air  density.  The  primary  energy  is  1  MeV.  Secondary  electron  produc¬ 
tion  has  been  restricted  to  less  than  20  keV. 


constituents  for  high  levels  of  energy  deposition.  Based  on  the  results  in 


Figure  S  which  apply  to  the  extremes  in  composition,  we  can  state  that 
conpositional  changes  one  might  generate  in  air  chemistry  modeling  will  not 
seriously  affect  that  part  of  the  secondary  flux  spectrum  capable  of  producing 
ionization  and  excitation.  It  should  be  noted  that  the  results  below  ~  10  eV 
are  sensitive  to  the  amount  of  plasma  present  through  the  derivative  term  in 
Equation  (2).  Furthermore,  it  is  worth  noting  that  the  most  pronounced 
difference  between  the  two  spectra  which  occurs  at  ~  2.5  eV  is  caused  by 
differences  in  vibrational  excitation.  For  Nj,  this  produces  strong  energy 
loss.  We  shall  limit  our  presentation  of  secondary  flux  results  to  those 
shown  in  Figure  5.  The  shape  of  the  flux  in  air  for  other  primary  energies  is 
similar  to  that  shown  for  Nj  over  the  given  energy  range. 

Before  proceeding  to  discuss  how  the  deposited  energy  becomes  distributed 
among  the  various  states,  we  wish  to  add  a  few  comments  to  those  given  in 
Section  7  on  Awip  •  Its  explicit  form  is 


A  W 


dB/dx 


lp  W:>b 


(13) 


where  dE/dx  is  the  stopping  power  and  Pion  is  the  sum  of  Pg  and  Pp.  dE/dx 
here  refers  to  only  the  deposited  energy  explicitly  treated  in  Equation  (2) . 
Thus,  it  is  less  than  obtained  by  the  Bethe  formula  if  the  maximum  secondary 
electron  energy  considered  is  smaller  than  its  maximum  allowed  value,  m  that 
case,  we  have  ussd  the  following  expression  for  dE/dx: 


dE/dx 


“*[k  M** 


W 


E  0,(E.  ,E  )dE 
8  l  O  S 


8 


(14) 


where  E  refers  to  the  largest  secondary  energy  treated,  n^  applies  to  N2  and 

02,  and  H^  refers  to  the  potential  energy  of  both  excited  and  Ion  states.  He 

wish  to  stress  that  the  needed  stopping  power  Is  the  one  just  discussed  since 

It  refers  to  the  actual  energy  deposition  being  modeled. 

Table  1  gives  A  w  for  the  four  primary  energies  considered.  For  the 

first  3  cases,  dE/dx  refers  to  total  energy  deposition  since  the  calculations 

were  performed  over  the  full  secondary  electron  energy  range.  Dissociative 

ionization  as  well  as  direct  ionization  have  been  included  in  Pn  and  P_ .  He 

P  3 

observe  some  variation  in  A  H.  but  in  all  cases  its  value  is  close  to  the 

ip 

accepted  value  of  ~  35  eV/ion  pair.  There  are  two  possible  sources  of  error 

in  the  calculated  values  of  A  H.  .  One  is  inaccuracies  in  the  calculated 

ip 

flux  vne(B).  The  other  is  in  the  applied  set  of  cross  sections  at  primary 
energies  needed  to  specify  dE/dx  and  Pp.  If  there  are  inconsistiencies 
between  this  set  and  its  counterpart  at  lower  energies  used  in  Equation  (2), 
they  will  affect  the  resulting  value  of  A  H  .  He  do  observe  some  error  in 
the  calculated  secondary  flux  based  on  the  applied  energy  conservation  test. 
The  error  is  not  serious,  however,  since  the  test  generally  leads  to  errors  of 
less  than  10%. 

18 

He  note  that  at  1  MeV,  P_  and  P_  are  comparable  whereas  Hirsh  et  al.  , 

P  3 

for  example,  report  relatively  larger  contributions  from  Pg  in  this  energy 

range.  They  observe  P  to  be  ~  2  P  .  He  would  also  obtain  approximately 

*  P 

this  result  had  we  allowed  for  the  energy  degradation  of  the  high  energy 
secondaries  extending  up  to  500  keV.  Both  Pg  and  dE/dx  would  increase 
accordingly  leaving  AH^  unchanged. 

He  now  address  the  question  of  what  the  distribution  of  energy  is  among 
various  states  of  interest.  What  will  be  shown  is  the  production  rate  of 
various  states  relative  to  the  total  ion  production  rate  P^on*  Table  2  gives 
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Table  X  —  Calculated  AW^  Values  and  Terms  Used  to  Specify  Them 


Ep(keV) 

(«V/cm) 

VJb 

Vib 

E  V^b 

AW.  (aV/ion  pair) 

1 

1.3(5) 

2.4(3) 

1.2(3) 

3.6(3) 

36 

10 

2.4(4) 

3.6(2) 

3.4(2) 

7.1(2) 

34 

100 

4.7(3) 

5.7(1) 

7.2(1) 

1.3(2) 

36 

1000 

1.4(3) 

2.3(1) 

1.9(1) 

4.1(1) 

34 

Table  2  —  Production  Given  Relative  to  Total  Ionization  Production. 

The  primary  energy  is  10  keV.  Similar  values  are  obtained  at  the  other  primary  energies 

considered  in  this  study. 


this  information  broken  down  into  direct  production  by  the  primary  electrons 
and  production  by  secondary  electrons.  The  results  apply  to  a  primary  energy 
of  10  keV.  We  have  generated  these  results  at  the  other  energies  as  well  and 
obtain  similar  values.  We  observe  similar  contributions  to  the  total 
production  by  primaries  and  secondaries  for  the  dipole  allowed  transitions 
producing  the  given  states.  For  triplets  states,  such  as  N2(A3£)  ,  the 
excitation  cross  section  falls  too  fast  with  increasing  energy  for  primary 
excitation  to  be  important.  The  total  values  for  ionization  and  dissociation 
in  Table  2  may  be  compared  with  similar  ones  generated  by  F.  Gilmore,  as 
quoted  by  Johnston  ,  using  a  continuous  slowing  down  scheme.  We  observe  no 
serious  differences  between  the  two  sets  of  values. 

9.  SUMMARY 

We  have  presented  a  computational  model  for  solving  an  integral  equation 
giving  the  secondary  electron  distribution  function  in  the  local 
approximation.  Qiergy  loss  is  treated  as  discrete  which  has  led  us  to  apply  a 
detailed  set  of  electron  impact  cross  sections.  Our  application  of  the  model 
has  been  to  energy  deposition  in  air  caused  by  the  passage  of  energetic 
primary  electrons.  Results  included  the  spectrum  of  secondary  electrons  under 
steady  state  conditions,  resulting  eV/ion  pair  values,  and  the  distribution  of 
energy  among  several  states.  The  calculated  eV/ion  pair  values  are  close  to 
the  universally  applied  value  of  ~  35.  The  distribution  information  presented 
gives  one  the  production  of  a  given  state  or  species  by  either  type  of 
exciting  electron- primary  or  secondary. 
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